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Abstract 

A recently introduced scheme for the renormalization of the Hartree-Fock-Bogoliubov equations 
in the case of zero-range pairing interaction is extended to the relativistic Hartree-Bogoliubov 
model. A density-dependent strength parameter of the zero-range pairing is adjusted in such a way 
that the renormalization procedure reproduces the empirical 1 S$ pairing gap in isospin-symmetric 
nuclear matter. The model is applied to the calculation of ground-state pairing properties of finite 
spherical nuclei. 
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I. INTRODUCTION 


The theoretical framework of self-consistent mean-field models enables a description of the 
nuclear many-body problem in terms of universal energy density functionals. By employing 
global effective interactions, adjusted to empirical properties of symmetric and asymmetric 
nuclear matter, and to bulk properties of spherical nuclei, the current generation of self- 
consistent mean-field models has achieved a high level of accuracy in the description of 
ground states and properties of excited states in arbitrarily heavy nuclei, exotic nuclei far 
from /3-stability, and in nuclear systems at the nucleon drip-lines |1]. 

The relativistic mean-field (RMF) models, in particular, are based on concepts of non- 
renormalizable effective relativistic field theories and density functional theory. They have 
been very successfully applied in studies of nuclear structure phenomena at and far from 
the valley of /3-stability. For a quantitative analysis of open-shell nuclei it is necessary 
to consider also pairing correlations. Pairing has often been taken into account in a very 
phenomenological way in the BCS model with the monopolc pairing force, adjusted to the 
experimental odd-even mass differences. This approach, however, presents only a poor 
approximation for nuclei far from stability. The physics of weakly-bound nuclei necessitates 
a unified and self-consistent treatment of mean-field and pairing correlations. This has led to 
the formulation and development of the relativistic Hartree-Bogoliubov (RHB) model, which 
represents a relativistic extension of the conventional Hartree-Fock-Bogoliubov framework, 
and provides a basis for a consistent microscopic description of ground-state properties of 
medium-heavy and heavy nuclei, low-energy excited states, small-amplitude vibrations, and 
reliable extrapolations toward the drip lines (2]. 

In most applications of the RHB model [2] the pairing part of the well known and very 
successful Gogny force []| has be employed in the particle-particle (pp) channel: 


V pp (l,2) = e 
1 = 1,2 


(( ri — r 2 )/Vi ) 2 (yy. + p.pv _ H .pr _ M .p«pT ) } (!) 

i=l,2 

with the set D1S 4] for the parameters //$, W % , Rj, Hi, and M % (i = 1,2). This force has been 
very carefully adjusted to the pairing properties of finite nuclei all over the periodic table. 
In particular, the basic advantage of the Gogny force is the finite range, which automatically 
guarantees a proper cut-off in momentum space. However, the resulting pairing field is non¬ 
local and the solution of the corresponding Dirac-Hartree-Bogoliubov integro-differential 
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equations can be time-consuming, especially in the case of deformed nuclei. Another pos¬ 
sibility is the use of a zero-range, possibly density-dependent, h-force in the pp-channel of 


0 . 


the RHB model |5l. This choice, however, introduces an additional cut-off parameter in 
energy and neither this parameter, nor the strength of the interaction, can be determined in 
a unique way. The effective range of the interaction is determined by the energy cut-off, and 
the strength parameter must be chosen accordingly in order to reproduce empirical pairing 
gaps. 


In a series of recent papers p, |7|, |8j A. Bulgac and Y. Yu have introduced a simple 
scheme for the renormalization of the Hartree-Fock-Bogoliubov equations in the case of 
zero-range pairing interaction. The scheme is equivalent to a simple energy cut-off with a 
position dependent coupling constant. In this work we use the prescription of Refs. 0, [?] 
to implement a regularization scheme for the relativistic Hartree-Bogoliubov equations with 
zero-range pairing. We analyze the resulting 1 Sq pairing gap in isospin-symmetric nuclear 
matter and apply the RHB model to the calculation of ground-state pairing properties of 
finite spherical nuclei. 

In Sec. d we present an outline of the RHB model and introduce the renormalization 
scheme for the case of zero-range pairing. The model is applied in Sec. IIIII to pairing in 
isospin-symmetric nuclear matter. Ground-state pairing properties of Sn nuclei are analyzed 
in Sec. EU Sec. ED contains the summary and conclusions. 


II. RELATIVISTIC HARTREE-BOGOLIUBOV MODEL WITH ZERO-RANGE 
PAIRING 

A detailed review of the relativistic Hartree-Bogoliubov model can be found, for instance, 
in Ref. 2]. In this section we include those features which are essential for the discussion of 
the renormalization of the RHB equations. The model can be derived within the framework 
of covariant density functional theory. When pairing correlations are included, the energy 
functional depends not only on the density matrix p and the meson fields 0 m , but in addition 
also on the anomalous density k 

^RHb\Pi k- (j)m\ E f{\[ J7 [p, (p rn ] T E pair [k] , (2) 
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where Ermf[p, 4 >] is the RMF energy density functional and the pairing energy E pair [k] is 
given by 


E pair [k] = -Tr [k*V pp k] . 


( 3 ) 


V pp denotes a general two-body pairing interaction. The equation of motion for the gener¬ 
alized density matrix 

/ 

p K 

-K* 1 — ft 


n = 


( 4 ) 


V 


reads 


id t U = [H(K),K] 


( 5 ) 


The generalized Hamiltonian EL is a functional derivative of the energy with respect to the 
generalized density 


'HRHB 


$Erhb 

<577 


^ h D — m — p, A 

v -A* -h* D + rn + p 


( 6 ) 


The self-consistent mean held ho is the Dirac Hamiltonian, and the pairing held reads 


A ab (r, r') = ^ V abcd( r > r>cd(r, r'), (7) 

Z c,d 

where a, b, c, d denote quantum numbers that specify the Dirac indices of the spinors, and 
Kifecd( r > r 0 are the matrix elements of a general two-body pairing interaction. 

Pairing effects in nuclei are restricted to an energy window of a few MeV around the 
Fermi level, and their scale is well separated from the scale of binding energies, which are in 
the range of several hundred to thousand MeV. There is no experimental evidence for any 
relativistic effect in the nuclear pairing held A. Therefore, pairing can be treated as a non- 
relativistic phenomenon, and a hybrid RHB model with a non-relativistic pairing interaction 
can be employed. For a general two-body interaction, the matrix elements of the relativistic 
pairing held read 


Aaipi,a2P2 r, E (aipi, a 2 p 2 \V pp \a 3 p 3 , a 4 p4)a ^a3P3,a4P4 j (8) 

^ a 3 p3,a 4 p4 

where the indices (p\ , p2, P:i, p4 — +> —) refer to the large and small components of the quasi- 
particle Dirac spinors. In most applications of the RHB model, only the large components 
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of the spinors £4(r) and 14 (r) have been included in the non-relativistic pairing tensor k in 
Eq. (THl . The resulting pairing held reads 

^ai+,a 2 + ~ ^ ^2 (°l + ’ a 2 + H PP \ a 3 + , «4+)a «a 3 + ,a 4 + • (9) 

a3+,a4+ 


The other components: A + _, A_ + , and A_are neglected, in accordance with the results 

that are obtained with a relativistic zero-range force [lOj]. 

The ground state of an open-shell nucleus is described by the solution of the relativistic 
Hartree-Bogoliubov equations 


/ 


hp — m — /i A 
-A* -h* D + m + n 



= El 



( 10 ) 


which correspond to the stationary limit of Eq. ©• 

The chemical potential /i is determined by the particle number subsidiary condition in 
order that the expectation value of the particle number operator in the ground state equals 
the number of nucleons. The column vectors denote the quasiparticle wave functions, and 
Eh are the quasiparticle energies. The RHB wave functions determine the hermitian single¬ 
particle density matrix 

Piv = (V*V T ) W , (11) 

and the antisymmetric anomalous density 

kiv = ( 1 v*u T )u >. ( 12 ) 

The calculated nuclear ground-state properties sensitively depend on the choice of the 
effective Lagrangian and pairing interaction. Over the years many parameter sets of the 
mean-held Lagrangian have been derived that provide a satisfactory description of nuclear 
properties along the /3-stability line. The most successful RMF effective interactions are 
purely phenomenological, with parameters adjusted to reproduce the nuclear matter equa¬ 
tion of state and a set of global properties of spherical closed-shell nuclei. This framework 
has recently been extended to include effective Lagrangians with explicit density-dependent 
meson-nucleon couplings. In a number of studies it has been shown that this class of global 
effective interactions provides an improved description of asymmetric nuclear matter, neu¬ 
tron matter and finite nuclei far from stability. In the present analysis of ground-state 
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properties of Sn isotopes the density-dependent effective interaction DD-ME1 {sj will be 
employed in the particle-hole (ph) channel of the RHB model. 

In the following we extend the regularization scheme of Bulgac and Yu [c], 3, Q] to the 
solution of the relativistic Hartree-Bogoliubov equations for a zero-range pairing interaction 


V pp (r, r') = g5 (r — r') . 


(13) 


In Refs. [fj 7] it has been shown that in this case the renormalized pairing held can be 
expressed as 

A(r) = -g eff ( r)re c (r) , (14) 

where K c (r) denotes the cut-off anomalous density 

E c 


K c(r) = Y V k( r )Uk(r) 

E k > 0 


(15) 


The cut-off energy E c dehnes the two corresponding momenta k c and l c 


k 2 (r) + m* 2 ( r) + R(r) — m = E c + g, 
\Jl 2 ( r) + m* 2 (r) + R(r ) — m = —E c + g . 


(16) 

(17) 


m*( r) = m+S( r) is the Dirac mass, and <S(r) and V (r) are, respectively, the scalar and vector 
single-nucleon potentials contained in the Dirac Hamiltonian ho- The chemical potential g 
determines the local Fermi momentum 


\Jk 2 f( r) + m* 2 ( r) + R(r) — m = g 
The effective, position-dependent coupling in Eq. (HI reads 

= - + Fi(r) + Fjfr) , 


(18) 


Self W a 


(19) 


with 


^i(r) = 


kf( r) . , k c (r) 

1 - Ar cth- v ; 


k c ( r) \Jk 2 f( r) + m* 2 (r) 

47r 2 

k c (r) \Jk 2 (r) +m* 2 (r) + m 


k c ( r) 


kfi 1 


8iT 2 


+ 


An 2 


kf(r) Ar cth 


2k 2 f {v) + ??7.* 2 (r) k c ( r) + Jk 2 ( r) + m* 2 ( r) 

J — In- 


87T 2 


m*( r) 


k c ( r) \Jk 2 {r) + m* 2 (v) 
kf(r)^k 2 (r) + m* 2 ( r) 

( 20 ) 
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Fo(r) = 


1 - k M Ar th- Ul) 


l c [r)^Jk 2 {r) + m* 2 ( r) 

47r 2 |_~ k( r ) 

l c (r)^Jl 2 c (r) + m* 2 (r) ^/fcj(r) +m 


M 1 


+ 


87T 2 47T 2 

2fc?(r) + m* 2 (r) Z c (r) + J/ 2 (r) + m* 2 (r) 

--- In--- 


fc/(r) Ar th 


l c ( r) \Jk 2 {r) + m* 2 ( 
kf(r)yjl 2 (r) + m* 2 ( 


(21) 


87t 2 m*(r) 

fq + F 2 is the relativistic generalization of the corresponding correction to the coupling 
constant g, as defined in Eq. (16) of Ref. jf|. 


III. PAIRING PROPERTIES OF SYMMETRIC NUCLEAR MATTER 

A zero-range pairing interaction leads to a particularly simple expression for the gap 
equation in symmetric nuclear matter 


1 

9eff 


47T 2 


elk 


k 2 


i 2 


( 22 ) 


+ A 2 


\Jk 2 + m * 2 — ^/k 2 + in * 2 

The momenta k c and l c are determined by the cut-off energy E c Eqs. mm, and the 
effective coupling g e jj is defined in Eq. (JTHJ). In the left panel of Fig. Q] we display the 
density dependence of the resulting pairing gap in nuclear matter (dashed curve). The 
single-particle spectrum has been calculated with the relativistic effective interaction DD- 
ME1 [9], and the coupling constant of the zero-range pairing interaction Eq. (11311 g = —330 
MeV fm 3 is typical for the values used by Bulgac and Yu in their analyses. The pairing 

t ap is shown in comparison to the gap calculated with the effective Gogny interaction D1S 
| (dots). The corresponding single-particle spectrum has been computed in the Hartree- 
Fock approximation for the Gogny interaction. The density dependence of the two gaps is 
completely different. The pairing gap of the renormalized zero-range interaction increases 
uniformly with density, whereas the gap of the Gogny interaction display the characteristic 
maximum of sa 2.5 MeV at low density p = 0.03 — 0.04 fm -3 (corresponding to a Fermi 
momentum of approximately 0.8 fm -1 ) and decreases at higher densities. The bell-shaped 
form of the pairing gap as a function of the density was, in fact, obtained already more 
than forty years ago jl2]. This density dependence is not characteristic only of the phe¬ 
nomenological finite-range interactions, but is also obtained when the gap is calculated with 
bare nucleon-nucleon potentials adjusted to the empirical nucleon-nucleon phase shifts and 
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deuteron properties (for a recent review see Ref. 


: G). 


The decrease of the gap at Fermi 


momenta kf > 0.8 fm 1 simply reflects the repulsive character of the nucleon-nucleon inter- 


innpl 

fit 


action at short distances Llll. Of course there is no repulsive component in the zero-range 
force with constant coupling Eq. m, and the corresponding pairing gap displays the un¬ 
physical uniform increase with density. We notice, however, that in the range of densities 
shown in Fig. |H i.e. up to nuclear matter saturation density, the values of the pairing gap 
of the renormalized zero-range interaction are comparable with those of the Gogny pairing 
gap. As will be shown in the next section, this means that the renormalization scheme for 
the zero-range interaction with constant coupling can be safely applied to the calculation of 
pairing correlations in finite nuclei, provided an appropriate choice is made for the strength 
parameter g. 

On the other hand, there is no particular reason why the strength parameter g of the zero- 
range pairing interaction should be a constant. In fact, in many applications to finite nuclei 
an explicit density dependence is introduced, and in this way pairing correlations partially 


G 


include finite-range effects. For instance, in one of the first applications IM Bertsch and 
Esbensen used a density-dependent contact interaction, together with a simple energy cut¬ 
off, in a description of pairing correlations in weakly bound neutron-rich nuclei. They also 
compared the corresponding pairing gap in symmetric nuclear matter with the result of 
a Hartree-Fock calculation using the Gogny interaction. In the present anaysis we have 
adjusted a density-dependent strength parameter g(p) of the zero-range pairing interaction 
Eq. (g3D , in such a way that the pairing gap of the renormalized zero-range interaction 
Eq. ({ 221 ) , reproduces the density dependence of the Gogny pairing gap. The resulting density 
dependence can be approximated by the following analytic expression 

1 


9(P) = 


(23) 


a 0 + Gyp 1 / 3 + a 2 p 2 / 3 

with ao = —0.064 fm -2 , a± = 0.447 fm -1 , and a 2 = —3.693. The resulting pairing gap, 
displayed in the left panel of Fig. (solid line), is in very good agreement with the one 


f.Q, 


calculated using the Gogny interaction. A very similar procedure was employed in Ref. 
where the density dependence of the “bare coupling constant” g(p) was adjusted to a specific 
formula for the pairing gap in low-density homogeneous neutron matter. 

In the right panel of Fig. Q]we display the effective couplings g e g- calculated using the 
constant g = —330 MeV fm 3 , and the density dependent coupling of Eq. (1231) . The density 
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dependence of the two effective couplings is completely different. In order to prevent an 
unphysical growth of the pairing gap with density, the density dependence of the pairing 
strength Eq. (EH ensures that the effective coupling becomes weaker with increasing nucleon 
density. A very strong effective coupling in the low-density region produces a peak in the 
corresponding pairing gap shown in the left panel. On the other hand, g e ff calculated using 
the constant coupling increases in absolute value with density, i.e. the resulting pairing 
gap increases uniformly with density. However, rather similar values for the two effective 
couplings g e jj are calculated in the region of densities characteristic for the bulk of finite 
nuclei. One should not, therefore, expect very different results for the pairing properties of 
finite nuclei calculated with the zero-range interaction with constant coupling, or with the 
density-dependent coupling of Eq. EH- In the next section we will show that this is really 
not true in weakly-bound nuclei far from stability. 

The renormalization prescription must, of course, lead to a pairing held which is inde¬ 
pendent of the cut-off energy E c , if the latter is chosen large enough. This is illustrated in 
Fig. El where we plot the pairing gap, calculated using the density-dependent coupling of 
Eq. EH ■ f° r a number of characteristic values E c in the interval between 5 MeV and 60 
MeV. The pairing gap shows a weak dependence on the cut-off energy only for the two lowest 
values of E c . When the cut-off is increased beyond 10 MeV, the corresponding pairing gaps 
cannot be distinguished. Thus already for E c > 10 MeV the pairing gap of the renormalized 
zero-range interaction in symmetric nuclear matter converges. This is in agreement with the 
results obtained in the analysis of the pairing gap in homogeneous neutron matter Q. 

IV. GROUND-STATE PAIRING PROPERTIES OF SPHERICAL NUCLEI 

In this section the renormalization scheme is tested in the calculation of ground-state 
pairing properties of Sn isotopes. The DD-ME1 mean-held Lagrangian is employed for 
the ph channel, and the zero-range interaction Eq. EH is used in the pp channel. The 
renormalization procedure described in the previous section is carried out for the zero-range 
interaction with constant pairing strength g = —330 MeV frn 3 , and and for the density- 
dependent coupling of Eq. EH- In the latter case the density dependence of the pairing 
strength has been adjusted to reproduce the Gogny D1S pairing gap in symmetric nuclear 
matter. In the following we denote by RCC the case of the renormalized constant coupling, 
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and by RDDC the results obtained with the renormalized density-dependent coupling. 

While in the symmetric nuclear matter the Fermi momentum is always real (see Eq. (ED), 
in the surface region of finite nuclei it becomes imaginary. In Ref. j^] it has been shown that 
also in this case the renormalized anomalous density is real. The effective coupling g e g is 
still given by Eq. ED, but 


FiW = 


k c ( r) \kf(r)\' 2 + m* 2 (: 

47T 2 

k c (r)yjk^(r) + m* 2 ( r) ^J—\kf(r)\ 2 + m 


kf( r) A k c (r) 

1 “ TT~\ Ar ctg TTT 

kc{r) k f ( r)J 


87T 2 


+ 


47T 2 


kf( r) Ar ctg 


(—21 A; y (r) | 2 + m* 2 (r)) fc c (r) + yjk 2 c (r) + m* 2 (r) 


87T 2 


m*{ r) 


k c (r)^/-\k f (r)\ 2 + m* 2 ( r) 
kf(j)\Jkl(r) + m* 2 (r) 

(24) 


and 


F 2 (r) = 0 . 


(25) 


However, if either k c or l c becomes imaginary, the corresponding terms in the effective 
coupling should be omitted. 

The rate of convergence of the renormalization scheme is illustrated in Fig. |3] where, for 
the nucleus 114 Sn, we display the average pairing gaps and the pairing energies as functions 
of the cut-off energy E c . The average gaps shown in the left panel, are defined as 


< A N >— 


^hnlj A nlj^nlj 


^2mlj V 


2 ’ 
nlj 


(26) 


where f 2 ^- are the occupation probabilities of the neutron states in the canonical basis. 
Both the pairing gaps and the pairing energies converge already for E c > 10 MeV. We also 
notice that, even though the renormalized constant coupling and the renormalized density- 
dependent coupling lead to very different pairing gaps in symmetric nuclear matter, in 114 Sn 
they produce similar average pairing gaps and virtually identical pairing energies. 

The corresponding pairing fields as functions of the radial coordinate, and g e jj Eq. ED 
as functions of the density, are plotted in Fig. 0 for a series of values of the energy cut-off. 
In both cases the calculation of the pairing field and g e g shows convergence for E c > 10 
MeV. While the renormalized constant coupling and the renormalized density-dependent 
coupling produce very similar average pairing gaps and pairing energies, the dependence 
of the corresponding pairing fields on the radial coordinate is rather different. The RCC 
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pairing field (upper left panel) is concentrated in the bulk of the nucleus, whereas the RDDC 
pairing held (lower left panel) exhibits a pronounced peak on the surface. This behavior 
reflects the difference between the effective couplings g e jj, already shown in the right panel 
of Fig. U for the case of symmetric nuclear matter. In the panels on the right of Fig. 0 
we plot the effective couplings g e ff( r (p)) as functions of the nucleon density in 114 Sn. The 
g e ff which corresponds to the density-dependent coupling of Eq. (1231) decreases steeply in 
the region of very low density, i.e., on the surface of the nucleus. Consequently, also the 
pairing held displays a peak in the surface region. In both the RCC and RDDC cases the 
pronounced discontinuity of the effective coupling g e jj at very low density corresponds to 
the transition from real to imaginary Fermi momentum kf. This is illustrated in Fig. 0 
where we plot the effective single-nucleon potential (left panel) and the correction to the 
coupling originating from the renormalization of the anomalous density (right panel). The 
effective single-nucleon potential is determined by the sum of the vector and scalar potentials 
Vcen(r) = S(r) + E(r). For real values of the Fermi momentum (the effective potential is 
below the chemical potential //) the correction to the coupling F\ (r) + F 2 (r) is calculated 
from Eqs. m and m, and for imaginary values of the Fermi momentum (the effective 
potential is above the chemical potential p) from Eqs. (12411 and (125(1 . In the region where 
the Fermi momentum changes from real to imaginary the correction Ei(r) + F 2 (r) displays 
a very sharp peak, which is reflected in the discontinuities of the effective couplings. 

The importance of possible surface effects is illustrated in Fig. El where we plot the 
calculated average pairing gaps and pairing energies for the chain of even-even Sn isotopes 
with 110 < A < 160. Although both the RCC and RDCC schemes lead to comparable 
values of the average pairing gaps for the entire isotopic chain, the pairing energies differ 
significantly for isotopes beyond the doubly closed-shell 132 Sn. For example, the pairing 
energy of 150 Sn calculated with renormalized density-dependent coupling (RDDC) is almost 
25 MeV larger than the one calculated with the renormalized constant coupling (RCC). The 
large increase in the pairing energy for the RDDC case is caused by the dominant role of the 
surface region for the very neutron-rich Sn isotopes, and because the effective coupling is 
especially strong at very low densities. In the panels on the left of Figs. [7] and |H1 we plot the 
self-consistent solutions for the cut-off anomalous densities Eq. m for the isotopes 114 Sn, 
124 Sn and 150 Sn, calculated using the RDDC and RCC effective couplings, respectively. The 
corresponding effective couplings g e jj are shown in the panels on the right of Figs. [3 and |H1 
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The anomalous densities for 114 Sn and 124 Sn are concentrated in the nuclear volume (r < 6 
fm), where the effective couplings g e jj have comparable values. Therefore, the corresponding 
pairing energies are similar for the RDDC and RCC cases. In 150 Sn, on the other hand, the 
anomalous densities extend to the region r > 8 fm, where the RDDC effective coupling 
becomes much stronger than the one calculated with the constant coupling (RCC). Hence, 
the pairing energy for 150 Sn, calculated using the renormalized density-dependent coupling 
is much larger than the one obtained with the renormalized constant coupling. 


V. CONCLUSIONS 


A simple renormalization scheme for the h 
range pairing has recently been introduced 


artree-Fock-Bogoliubov equations with zero- 
In the present work we have implemented 


u, 

this renormalization scheme for the relativistic Hartree-Bogoliubov equations with a zero- 
range pairing interaction. The procedure is equivalent to a simple energy cut-off with a 
position dependent coupling constant. We have verified that the resulting average pairing 
gaps and pairing energies do not depend on the cut-off energy E c . if the latter is chosen large 
enough. Convergence is achieved already for values E c > 10 MeV, both in nuclear matter 
and for finite nuclei. If the strength parameter of the zero-range pairing is a constant, 
the resulting pairing gap in symmetric nuclear matter displays an unphysical increase with 
density. We have therefore adjusted a density-dependent strength parameter of the zero- 
range pairing in such a way that the renormalization procedure reproduces in symmetric 
nuclear matter the pairing gap of the phenomenological Gogny interaction. In this sense the 
present study goes beyond the simple extension of the renormalization scheme of Ref. (J 
to the relativistic framework. However, the resulting effective coupling is too strong in 
the region of low density, and this leads to large pairing energies in open-shell nuclei with 
very diffuse surfaces, e.g. in neutron-rich Sn isotopes. One must therefore be careful when 
applying the renormalized HFB or RHB models with zero-range pairing to nuclei far from 
stability. Adjusting the strength parameter to the pairing gap in symmetric nuclear matter 
obviously does not provide enough information about the density dependence of the zero- 
range pairing to be used in very neutron-rich nuclei. 

ACKNOWLEDGMENTS 

This work has been supported in part by the Bundesministerium fur Bildung und 


12 




Forschung - project 06 MT 193, by the Alexander von Humboldt Stiftung, and by the 
Croatian Ministry of Science - project 0119250. 


[1] M. Bender, P.-H. Heenen, and P.-G. Reinhard, Rev. Mod. Phys. 75 (2003) 121. 

[2] D. Vretenar, A.V. Afanasjev, G.A. Lalazissis, and P. Ring, Phys. Rep. (2005). 

[3] J.F. Berger, M. Girod, and D. Gogny, Nucl. Phys. A 428, 23c (1984). 

[4] J. F. Berger, M. Girod, and D. Gogny, Commput. Phys. Commun. 63, 365 (1991). 

[5] J. Meng, Nucl. Phys. A 635, 3 (1998). 

[6] A. Bulgac, Phys. Rev. C 65, 051305 (2002). 

[7] A. Bulgac and Y. Yu, Phys. Rev. Lett. 88, 042504 (2002). 

[8] Y. Yu and A. Bulgac, Phys. Rev. Lett. 90, 222501 (2003). 

[9] T. Niksic, D. Vretenar, P. Finelli, and P. Ring, Phys. Rev. C 66, 024306 (2002). 

[10] M. Serra and P. Ring, Phys. Rev. C65, 064324 (2002). 

[11] M. Serra, A. Rummel, and P. Ring, Phys. Rev. C65, 014304 (2002). 

[12] V.J. Emery and A.M. Sessler, Phys. Rev. 119, 248 (1960). 

[13] D.J. Dean and M. Hjorth-Jensen, Rev. Mod. Phys. 75, 607 (2003); and references therein. 

[14] G.F. Bertsch and H. Esbensen, Ann. Phys. 209, 327 (1991). 

[15] Y. Yu and A. Bulgac, Phys. Rev. Lett. 90, 161101 (2003). 


13 




FIG. 1: Pairing gap in symmetric nuclear matter as a function of density for the zero-range 
interaction with constant coupling g = —330 MeV fm 3 (dashed) and the density-dependent coupling 
Eq. (f?3l) (solid). The corresponding density-dependent curves g e ff (Eq. (TlTflO are plotted in the 
panel on the right. The dots in the left panel denote the pairing gap calculated with the Gogny 
D1S interaction. 
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FIG. 2: Pairing gap in symmetric nuclear matter as a function of the Fermi momentum, calculated 
with the zero-range interaction and the density-dependent coupling Eq. (El, for a series of cut-off 
energies. 
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E c (MeV) E c (MeV) 

FIG. 3: Average neutron pairing gaps (left), and pairing energies (right) for 114 Sn, calculated 
with a zero-range interaction, as functions of the cut-off energy. The calculations are performed 
with the renormalization of the constant coupling g = —330 MeV fm 3 (squares), and with the 
renormalization of the density-dependent coupling Eq. (f?3l) (dots). 


16 














6 

r(fm) 


8 10 



0.12 


0.12 


FIG. 4: The pairing fields as functions of the radial coordinate (left), and the curves g e g(r(p)) 
(Eq. (IIHD ) as functions of the corresponding density in 114 Sn (right), for a series of energy cut¬ 
offs. The upper panels display results of the renormalization procedure for the zero-range force 
with constant coupling g = —330 MeV fm 3 (RCC), and the lower ones correspond to the density- 
dependent coupling Eq. ESI ) (RDDC). 
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FIG. 5: The effective single-nucleon potential in 114 Sn as a function of the radial coordinate (left), 
and the correction to the strength parameter of the zero-range effective force (right) as a function 
of the square of the Fermi momentum. // denotes the position of the chemical potential. 
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FIG. 6: Average neutron pairing gaps (left), and pairing energies (right) for the chain of even-even 
Sn isotopes with 110 < A < 160. The calculations are performed with the renormalization of 
the constant coupling g = —330 MeV fm 3 (dashed), and with the renormalization of the density- 
dependent coupling Eq. (12311 (solid). 
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FIG. 7: The anomalous densities in 114 Sn, 124 Sn, 150 Sn isotopes as functions of the radial coordinate 
(left), and the density-dependent curves g e g (Eq. (1191) 1 (right). The calculations are performed 
with the renormalization of the density-dependent coupling Eq. pill (RDDC). 
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FIG. 8: Same as in Fig. 0 but for the RCC case with the constant coupling g = —330 MeV fm 3 . 
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